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Abstract. We present spectroscopic observations of the runaway reddened star HD 34078 acquired during the last three years 
at Observatoire de Haute Provence and McDonald Observatory as well as other spectra obtained since 1990. The drift of the 
line of sight through the foreground cloud due to the large transverse velocity of HD 34078 allows us to probe the spatial 
distribution of CH, CH+, CN and DIBs carriers at scales ranging from about 1 AU up to 150 AU. In particular, time variations 
in the equivalent width of absorption lines are examined. A few past and recent high resolution observations of CH and CH+ 
absorption are used to search for line profile variations and to convert equivalent widths into column densities. 
The data set reveals a 20% increase in CH column density over the past 10 years with no corresponding variation in the column 
density of CH+ or in the strengths of the 5780 and 5797 A DIBs. CN observations indicate that its excitation temperature 
has significantly increased from < 3.1 K in 1993 to 3.6 ± 0. 17 K in 1998 while the CN column shows only a modest rise of 
~ 12 ± 6%. The data also strongly suggest the existence of weak correlated variations in CH and CH+ columns over periods 
of 6 - 12 months (or x 10 AU). 

These results are discussed in relation to CH+ production mechanisms. A dense newly intervening clump is considered in 
order to explain the long-term increase in the column density of CH, but such a scenario does not account for all observational 
constraints. Instead, the observations are best described by CH+ production in a photodissociation region, like that suggested 
for the Pleiades and IC 348. 
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1. Introduction 

While the existence of AU-scale structure is relatively 
well established in atomic gas (see e.g., pieter et al., 1976 
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iFaison & Goss, 2001 
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[Welty & F itzpatrick, 200]), the reality 
of such tiny fluctuations within molecular gas is still ques- 
tionable. Indeed, only minor species can be easily observed 
towards molecular clouds and their spatial distribution might 
not reflect that of H2 if "chemical structure" were present. 
H2CO, HCO+ and OH apparently display column density fluc- 
tuations reaching 5 to 15% along lines of sight separated by 
about 10 AU (Marscher ^tal., 19 93 Moore & Marscher, 1995 
ILiszt & Lucas, 2000j , while dust grains appear to be more 
smoothly distributed fThoraval et al., 1996 1. At larger scales 
(about 10,000 AU), Pan et al. (2001) find significant differences 
in CN, CH and CH+ absorption lines. 

If marked enough, small scale structure within molecular 
gas might notably affect its time evolution, fragmentation and 
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* Based on observations made at Observatoire de Haute Provence 
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then utimately, star formation. It is therefore important to char- 
acterise the properties of such media, identify the parameters 
displaying fluctuations and quantify their amplitude and scale- 
length. To this aim, we have selected a bright 09.5 runaway 
star, HD 34078 (AE Aur), seen through a translucent cloud 
with E(B -V) = 0.52 (.Diplas & Savage, 1994| l. Thanks to its 
large transverse velocity, about 100 km s~^ for an assumed dis- 
tance of 530 pc, the comparison of spectra taken at one year in- 
tervals provides a measurement of column density variations 
at a scale of 17 AU (assuming a cloud distance of 400 pc: 
IBrown et al., 1995l l. To complement ongoing FUSE observa- 
tions designed to investigate in a direct way small scale varia- 
tions of the H2 column density (and then study density structure 
within molecular gas), we have undertaken repeated ground- 
based observations of CH, CH+and CN absorption lines and 
diffuse interstellar bands (DIBs) towards HD 34078. 

These optical spectra complement the FUSE UV data in 
several important respects: 

- since CH is seen to correlate well with H2 over large 
scales (e.g., |Federman, 1982t , it should be a good indica- 
tor of cloud structure. On the other hand, the abundance of 
CH+ is largely independent of H2 column density, A^(H2); the 
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large abundance of this species is still poorly understood and 
hot pockets of gas (e.g., shocks or vortices as proposed by 
Falgarone & Puget, 1995) could be the formation site. Then, 
CH+ observations may reveal the structure of these presumably 
very small regions where energy is actively dissipated. Finally, 
CN depends nonUnearly on H2 and is thus an indicator of gas 
density ( |Federman et al., 1994t . 

- given the brightness of HD 34078, high S/N spectra can 
be obtained easily to search for very small fluctuations which 
would be difficult to detect in FUSE spectra, 

- spectral data in the visible range can be acquired much 
more easily than in the far UV, and therefore, it is possible 
to control the time (i.e., spatial) sampling and to accumulate 
a larger number of spectra. Moreover, several visible observa- 
tions of HD 34078 have been performed in the past 10 years, 
thus giving access to large scale variations. A broad range of 
scales can then be explored through visible observations. 

- high resolution observations providing information on the 
velocity distribution and its possible variations are feasible in 
the visible, but not in the FUSE range. 

In this paper, we present observations made specifically for 
this project together with a comparison with older data. The 
whole set of spectra is used to search for time variable absorp- 
tion and infer the structure implied for CH, CH+, CN and dif- 
fuse bands carriers. Sect. 2 presents the observations and the 
methods used to analyse the data. In Sect. 3, we give the results 
obtained for CH, CH+, CN and some selected DIBs. Since 
most observations have been obtained at low resolution, we 
mainly consider equivalent width variations. High resolution 
observations of CH and CH+ are then used to infer velocity 
distributions and translate variations in equivalent width into 
column density variations for these two species. The implica- 
tion of our results for the existence of AU-scale structure in the 
distribution of these tracers is discussed in Sect. 4. 

The analysis of Cycle 1 and Cycle 2 FUSE data wiU be pre- 
sented in a separate paper (Boisse et al., 2002; see also Boisse 
et al., 2001 and Le Petit et al., 2001 for a preliminary report). 
Our ultimate goal is to coiTelate the variations seen for H2 and 
those of other tracers like CH, CH+ and CN in order to improve 
the description of the structure. 

2. The Data 

2. 1. The Available Spectra 

Spectra acquired at various observatories have been used for 
this study. In the following, these are named according to the 
observer and year of data taking. We obtained specifically for 
this project a series of nine spectra at the Observatoire de Haute 
Provence (OHP) (R99-02). In addition, high resolution spectra 
were recently recorded at McDonald Observatory (F02). The 
motivation was twofold. First, as we are interested in column 
density variations, we need to translate the observed changes 
in equivalent width, W\, into changes in N (i.e., determine N 
and dN/dW\). Since the lines are not optically thin, the con- 
version requires a good knowledge of the velocity distribution. 
This point is especially important for CH A4300 because this 
feature is in fact a blend of two transitions and the analysis 
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Fig. 1. The impact point of the line of sight towards HD 34078 
onto the foreground cloud versus time. The motions of the 
earth, the sun and HD 34078 have been taken into account. 
Markers specify the position of the star for each epoch at which 
spectra were obtained (cf. Table IaTH . Note that the X and Y 
scales are not identical. The insert focuses on the recent obser- 
vations made for this project (R99a-R02, F02) that probe the 
smallest scales. 



can be notably affected by the splitting if the b parameter is 
less than 2 km s~^ ( |Lien, 1984) . Second, only high resolution 
observations allow us to determine whether or not line shapes 
have changed, and they possibly reveal profile variations not 
associated with equivalent width changes. 

HD 34078 was observed several times since 1991 by 
one of us (S. Federman at McDonald Observatory: F93a, 
F93b). Spectra of HD 34078 were also acquired by M. Allen 
(KPNO: A91), P Jenniskens and FX. Desert (at OHP: J91), 
J. Krelowski, G. A. Galazutdinov & F. A. Musaev (at Terskol 
Observatory: K97a, K97b) and G. Herbig (at the Keck tele- 
scope: H98) who kindly provided their data, so that we could 
analyze them in the same way as our own and search for varia- 
tions by direct comparison of the spectra. Note that in his study, 
Herbig (1999) searched for changes in line position with re- 
spect to the spectrum obtained by Adams (1949) and thus in- 
vestigated the velocity structure of the foreground gas; no sig- 
nificant variation was observed. 

Table lA. II summarises all these observations and gives the 
date, resolution and signal-to-noise ratio. The transitions (from 
CH, CH+, CN or DIBs) contained in the observed wavelength 
range are indicated. These data altogether provide a sort of ID 
cut through the cloud, the sampling of which is illustrated by 
Fig. n Below, we give some details concerning the observa- 
tions of HD 34078 used in this paper and not previously pub- 
hshed. 
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Fig. 2. The CH+ lines at 3957 and 4232 A (a,b) and CH line at 
4300 A (c) in the OHP spectrum recorded in December 1999. 
The raw spectrum has just been divided by a constant to scale 
the mean intensity to 1. Note the broad Fe II stellar feature 
blended with the 4232 A line. The S/N is 230 per 0.03 A pixel 
around 4300 A. 

2.2. OHP Observations 

The spectra were obtained using the fiber-fed echelle spec- 
trograph Elodie ( Bara nne et al., 199'6t mounted on the 1 .93 m 
telescope (R99-02). The resolution is about 32000. Initially, a 
spectrum was taken each month in order to probe very small 
scales. Since this series did not reveal significant variations, 
only two spectra were taken during the next winters. The range 
covered is 3906 A - 681 1 A and includes lines from CH (4300 
A), CH+ (3957 and 4232 A) and DIBs. For all spectra, the 
integration time is about 1 hour, split in 2 to 3 successive ex- 
posures of 30 or 20 minutes each in order to check the stability 
of the instrument and robustness of our analysis on indepen- 
dent spectra. The corresponding S/N ratio is typically 150 per 
resolution element. Some examples of CH and CH+ lines de- 
tected in OHP spectra are shown in Fig. |2] The spectra have 
been extracted and wavelength calibrated using the automatic 
on-line data reduction program attached to the instrument (see 
IBaranne et al., 1996| for details); the wavelength calibration is 
based on Thorium lamp spectra. 

2.3. McDonald Observations 

Three data sets were acquired since the last published results 
(Federman et al., 1994) by one of us (S.R.F.). Two of the new 
sets utilized the 6 foot camera and the 2dcoude spectrograph 
on the 2.7 m telescope. A standard observing strategy, with 
bias frames, flat fields, and a Th-Ar comparison spectrum, was 
employed (see Knauth et al., 2001 for details of similar mea- 
surements). Individual orders containing CH+ A4232 and CH 
A4300 were imaged with the 6 foot camera in 1993. The nomi- 
nal resolution was 1.5 km s~^. In 2002, spectra of CN, CH and 
CH+ were acquired at a single setting of the 2dcoude spec- 
trograph. Here, the resolution was about 1.9 km s^^. Also in 
1993, the Sandiford echelle spectrograph on the 2.1 m tele- 
scope (McCarthy et al., 1993) was used for measurements on 



various species in a similar fashion. One setting, centered on 
3990 A, provided data on CN, CH AA3878,3886,3890 and 
CH+ A3957. The second setup at 4300 A measured absorption 
from CH+ A4232 and CH A4300. The widths of Th-Ar fines 
indicated a resolution of ^ 7.5 km s^^ for these spectra. All 
data were reduced and analyzed with the IRAF package. 



2.4. Analysis of the Data 

In most of the available HD 34078 spectra, the profiles are not 
resolved. We shall then look primarily for variations in W\. 

We have been particularly careful in placing the continuum 
in a consistent manner for all epochs and also in estimating 
properly the uncertainty in W\. For some species, several lines 
are observed, e.g., at 3957 A and 4232 A for CH+. Similarly, 
a few observations of the CH fines at 3878, 3886 and 3890 
A are available for comparison with the 4300 A results. This 
redundancy is useful as a means to assess the robustness of our 
analysis and validate the error estimates. 

We have developed an automatic procedure based on the 
MIDAS package alice in order to fit the continuum in a user- 
independent way and to measure W\ values. Several "clean" 
windows (generally two, located on both sides of the line) are 
selected; a polynomial fit with a degree ranging from 1 to 3 and 
depending on the shape of the spectrum is then performed. This 
provides a good estimate of the continuum in the interval cov- 
ered by the absorption feature (see Fig.|2}. One cr uncertainties 
in W\ can be estimated as the quadratic sum of e„ and Ec, the 
pixel-to-pixel noise and the uncertainty in the continuum level 
in the normafised spectrum respectively. The latter are given by 



— (Tp X SX X V N , and 
e, = AA,,;„,. X I 1 - F I , 



(1) 
(2) 



where AA/i„e is the total width of the absorption feature, SX is 
the pixel size and N the number of pixels covered by the line 
{AXiine = NSX). The average (Eq.|3 and the rm.s. (Eq. Q 
of the normalised flux, F, are computed along the windows 
adjacent to the line. The robustness of the fitting procedure and 
the estimates for uncertainty have also been checked against 
synthetic spectra including absorption lines with known Wx, a 
signal-to-noise ratio comparable to that in our data, and various 
continuum shapes. 

Some transitions present difficulties. CH+ A4232 is blended 
with a shallow stellar feature (from Fe II, see Fig.|2j3) and the 
above procedure cannnot be applied. Similar problems are en- 
countered with the R(0) CN line at 3876 A. In such cases, 
the W\ is measured using the command integrateAine from 
MIDAS (an order 2-3 is chosen for the polynomial used to fit 
the spectrum immediately adjacent to the line) and the uncer- 
tainties are determined from successive measurements. 

The high resolution profiles have been fitted using VPFIT 
( ICarswell et al., 1987' l. No attempt was made to look for veloc- 
ity shifts (the question addressed by Herbig, 1999). 
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3. Results 

We now present the results obtained for CH, CH+, CN and 60 
some selected DIBs. For CH and CH+, we first discuss the ho- 
mogeneous set of data obtained after November 1999 which 
probe time scales ranging from one month to 2.5 years (cor- 
responding to 1.5 - 40 AU). The whole set of available obser- 
vations made in the past 1 1 years are presented next. Finally, g 
we discuss the high resolution data and their implications for 
velocity structure and column density variations. 



50 



3. 1. structure in CH 

3.1 .1 . Small scale variations 

Fig.|3t displays the evolution of the equivalent width as a func- 
tion of time for the 4300 A line. The variation since 1999 is 
fitted with a linear form W\ = W\,o + a{t ~ tmcd) where 
tmed is the median of observing times (December 2000). The 
data points are weighted according to signal-to-noise ratio in 
all fits performed in this paper. The contours of the reduced 
X^(Wo, a) are plotted in Fig. 0^ (dashed lines) for the linear 
case. A constant value is rejected only at the la level. The best« 
fit, Wq = 56.6 mA and a = 0.98 mA yr~^ corresponds to - 
AWx ^ 2.2 mA between November 1999 and February 2002. : 
It is overplotted in Fig.|3t (dotted line). Large variations at this 
scale are excluded since —1.8 < AW\ < 6.3 mA at the So- 
level. 

The linear fit is not fully satisfactory since residuals corre- 
spond to 2(T. This suggests that small additional eiTatic varia- 
tions in W\ may be present. The observed amplitude is 1.50 
mA (rms) about the linear fit (and 1 .56 mA about a constant 
value; the maximum difference is 4.8 mA). A Kolmogorov- 
Smirnov test (Miller, 1956) shows that these values are too 
small compared to observational errors to indicate with cer- 
tainty that we are detecting fluctuations in Wa(CH). However, 
the reality of the latter is supported by the fact that similar vari- 
ations are observed for CH+ lines (see below). 

3.1 .2. Large scale variations ; 

Results for the full set of observations are given in Table IX2I 
and shown also in Fig. [3^. A significant increase is appar- 
ent. Fitting again with a linear variation (with imod = April 
1996) gives a total increase of AW\ = 8 mA or 16% between 
1992 and November 2001, consistent with the slope derived 
from the small scale data alone (see Fig. |4j. The 3a limit is 
4.2 < AWx < 13 mA. However, the fit is even less satisfac- 
tory than for the recent data alone since the minimum is 
larger than 3 a. The assumption of a linear increase is some- 
what arbitrary and visual inspection of Fig. |3] suggests that a 
step function might also provide an acceptable fit. The mini- 
mum is not significantly lowered, ~ 2.6a. The date at which 
the step occurs is well constrained by K97 and H98 measure- 
ments. The amplitude (AVFa ^ 7 mA) is similar to the total 
increase infeiTed from the linear fit. To assess the reality of an 
increase in Wa(4300) regardless of any assumption on its form, 
we performed a Pearson test to investigate the correlation be- 
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Fig. 3. The equivalent width of CH and CH+ lines versus 
time. Top panel: CH A4300; middle panel: CH+ A4232. bottom 
panel: CH+ A3957. Solid markers show recent closely-spaced 
OHP observations. A linear fit is performed on the whole data 
set (solid lines). For recent OHP/F02 data on CH A4300, a sep- 
arate linear fit is given (dotted line on top panel). A fit by a step 
function is also shown for the CH data, with eiTors for the upper 
and lower values indicated by the small vertical bars. Note the 
similarity between CH and CH+ variations from October 2000 
to March 2002 (points labelled A, B, C and D), suggesting the 
reality of correlated erratic fluctuations. Such fluctuations are 
also suggested by the McDonald data of F02. 
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Fig. 4. Linear fit of Wx{t) for CH A4300 (a) and CH+ A4232 
(b) (Wo = WA(tmcd))- Solid line: linear fit to the whole set 
of data (with imod = April 1996). Dashed line: linear fit to the 
recent data (tmcd = December 2000). The best fits are shown in 
Fig- 13 Contours represent the levels 1, 2 and 3 of the reduced 
(y^ )• 
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file was fitted using VPFIT, taking into account the doublet 
structure of the CH ground level (Lien, 1984). Both levels are 
assumed to have the same column density. Three components 
are needed: two narrow ones reproduce the slightly asymmetric 
core of the absorption line and a much broader and fainter one 
accounts for the wings which are present on both sides of the 
core. The resulting fit and the three components are shown in 
Fig. |5] and the parameters of the fit are given in Table [2 There 
is no doubt about the reality of the broad component since it is 
clearly seen in all spectra with the appropriate S/N ratio (A91, 
F93a, F02) and for both the 3886 and 4300 A lines (H98). It 
is noteworthy that, although barely noticeable, it contains no 
less than 28% of the total CH column density. As a consistency 
check, we verified that W\ values measured for all CH lines 
detected in Herbig's spectrum are consistent with the velocity 
distribution derived from the F02 spectrum. We note, however, 
that different 3-component solutions for the latter are also ac- 
ceptable (in particular, one in which the strongest narrow com- 
ponent is the red one). Nevertheless, since the optical depth 
remains moderate, this degeneracy does not result in larger un- 
certainties on iV(CH). Moreover, very recent observations of 
the CO(2 - 1) emission line at 1.3 mm towards HD 34078 per- 
formed at the 30 m IRAM telescope give a velocity profile in 
excellent agreement with that shown in Fig.|5](E. Roueff and 
M. Gerin, private communication). 



4300 



4300.5 



Fig. 5. High resolution observation (F02) of the CH line, solid 
line: fit with VPFIT. dotted lines: the three components of the 
fit. 



tween W\ and time. An uncorrected evolution of CH A4300 
with time is rejected at the 4 a level. 

Other CH lines at 3886 and 3890 A have been observed 
in 1990, 1993 and 1998. Their variation should be consistent 
with the increase seen in CH A4300. Given the weakness of 
these lines and the S/N attained, we find that even the strongest 
feature at 3886 A does not bring useful additional constraints. 

3.1.3. High resolution observations 

We first discuss the highest signal-to-noise ratio observation of 
CH A4300, performed in 2002 at McDonald (F02). The pro- 




4300 4300.5 

Wavelength (A) 

Fig. 6. Variation of the CH profile. Top panel: F93a and F02 
spectra are overplotted as thick and solid line respectively. 
Bottom panel: the optical depth of F93a is multiplied by 1.34 
and the resulting spectrum is subtracted from F02. The residual 
(solid line) is compared to the la noise of F93a (dotted line). 



If we now fit independently the F93a profile, a slightly dif- 
ferent solution is obtained. The parameters are less accurate 
due to the lower signal-to-noise ratio — see Table [2 The real- 
ity of profile variations is difficult to assess from comparison 
of the F93a and F02 parameters because the decomposition is 
not unique; we give in Table^solutions with the lowest but 
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Table 1. Relative velocity, Doppler parameter and column den- 
sities derived by fitting the high resolution CH and CH+ pro- 
files with VPFIT. Three components are assumed; results ob- 
tained with one single component are also given. Solutions with 
the lowest are listed but relatively different solutions cannot 
be rejected. The doublet structure of the CH ground level has 
been taken into account: the column density given in this table 
is the sum for the two levels. The reduced is always ^ 1. 



relatively different solutions cannot be rejected. We then per- 
form a more direct comparison of profiles and find that multi- 
pUcation of the F93a optical depth by 1.34 provides an accept- 
able fit to the F02 profile (Fig.|6j the slightly different spectral 
resolutions of the two spectra have been taken into account). 
Therefore, we conclude that the velocity distribution has re- 
mained roughly constant between 1993 and 2002 and that the 
observed line variation is essentially due to an increase by the 
same factor of the column density for each of the three compo- 
nents. 

For the velocity profile obtained above, the ratio 
(SN/Ny(6Wx/Wx) needed to translate equivalent width vari- 
ations into N variations increases from 1.37 to 1.50 between 
Wx = 48 and W\ — 60 mA. For a single component pro- 
file, these ratios equal 1 .36 and 1 .49 respectively, showing that 
in our optical depth regime, (6N/N)/(SW\/W\) depends very 
weakly on the exact model chosen to fit the profiles. This im- 



plies an increase of AN/N — 22% in either the step function 
fit or in the linear fit (between 1992 and 2001). In this latter 
case, bounds of 12% < AN/N < 38% ai-e obtained at the So- 
level. The higher (34%) increase in A^(CH) between F93a and 
F02 corresponds to a 12 mA increase in W\. 

3.2. Structure in CI-I+ 

CH+ ti-ansitions are observed at 3957 A and 4232 A in the OHP 
data; some additional lines are detected at shorter wavelengths 
in G. Herbig's spectrum. The absorption line at 3957 A is free 
of blending and the continuum can be accurately determined. 
Although the 4232 A line is more difficult to measure, the S/N 
ratio in this region is higher and thus comparable accuracy is 
obtained for the two lines. 

3.2.1 . Small scale variations 

Both lines appear to be roughly constant in equivalent width 
(Fig.|3j), c). Erratic variations still seem to be present, although 
the larger errors yield an acceptable value of '^^^ 1 for the best 
fit. In particular, the pattern characterising the time variations 
is similar for both CH+ lines and for the CH A4300 transition. 
This is especially clear for the last four OHP measurements, 
marked with letters A to D in Fig. |3] We also note that the 
F02 measurement, which has been made in completely inde- 
pendent conditions, apparently displays this coherent CH/CH"*" 
behavior Finally, it appears in Fig.|3t and Fig.lSj) that closely- 
spaced measurements (made during winter 1999-2000) display 
relatively smaller scatter, indicating a possible cut-off in the 
structure at scales below 1-2 AU. 

Such a coherent variation could be an artefact due to in- 
strumental effects (e.g., if changes in the background were not 
taken into account properly). We searched for any dependence 
on the observing conditions (e.g., high or low sky background) 
and found none. We also selected a few narrow stellar lines 
near 4300 A for which the continuum level is well defined and 
measured their equivalent width; no significant variations of 
the kind observed for CH or CH+ are observed. Finally, we 
checked that the depth of the stellar H7 absorption is stable as 
it should be in the absence of uncontroled variations of the zero 
level. 

We therefore strongly suspect that stochastic correlated 
fluctuations of 7V(CH) and A^(CH+) are present, although more 
accurate measurements (especially for CH+) would be useful 
to establish more firmly their reality and to characterise their 
properties. Among measurements performed after November 
1999, the r.m.s. fluctuation of W^a(4300) is 1.5 mA while that 
of PVa(4232) is 1.8 mA (the corresponding fluctuations in N 
are discussed below). 

3.2.2. Large scale variations 

Fig. shows the result of a linear fit to the whole set of 
data for CH+ A4232. The minimum is again close to 1. 
Although the best fit is a slow decrease, dW\/dt — —0.4 mA 
yr^^, a constant value, W\ ~ 40 — 43 mA, cannot be re- 
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Fig. 7. High resolution observation (F02) of the CH+ A4232 
line {dotted lines: the three components of the fit and the stellar 
line; solid line: resulting fit). 

jected at the 1.5 a level. We obtain the following 3(t bounds: 
-1.2 < dWx/dt < 0.45 mA yr"^ Clearly, the variation ob- 
served for CH+ is quite different from that seen for CH. 

Consistent results are obtained from CH+ A3957. The ratio 
of the two total CH+ equivalent widths, Wx{'i951)IWx{^232), 
is found to be constant at a value of 0.65 ±0.05, indicating that 
these lines are nearly optically thin (ratio of 0.60). The best 
linear fit is also decreasing, although a constant value for both 
Wa(3957) or W^x(4232) cannot be ruled out. 

3.2.3. High resolution observations 

CH+ A4232 has been observed twice at high resolution (F93a 
and F02). To fit these profiles, we first remove the blended stel- 
lar line using a specific Gaussian profile. The CH+ line in F02 
is asymmetric, as was the CH A4300 one, but now it is the blue 
wing which is more extended. A good fit is obtained with 3 
components; parameters are summarised in Tableland shown 
in Fig.El 

The F93a and F02 profiles are almost identical. The same 
components fit both observations with a residual less than 
1.5 and 1.0 cr, respectively. If the F93a spectrum is analysed 
separately, a good fit is obtained with two narrow components 
only, but this is just due to the lower signal-to-noise ratio. The 
third broad component is also apparent in the H98 spectrum 
for both the CH+ fines at 3957 and 4232 A. 

In summary, the CH"*" line profiles have been stable be- 
tween 1990 and 2002. For the model derived from the F02 pro- 
file and the CH+ A4232 data, one gets 5N/N =1.3* 5Wx/Wx 
in the range 5Wx = —10 to 10 mA. We then derive an up- 
per limit for the annual variation of the column density from 
1990 to 2002 of -3.9 < 5N/N < 1.45% per year at the 
3(7 level. Tighter constraints can be obtained by combining 
CH+ A3957 and CH+ A4232 measurements. For all epochs, 
we derive iV(CH"'") and (t{N) from each transition; estimates 
drawn from CH+ A4232 appear to be, in average, slightly 
lower than those inferred from CH+ A3957, presumably be- 



T 




LSR Velocity (km/s) 

Fig. 8. Comparison of optical depth profiles for high resolution 
observations (F02) of CH A4300 (solid line) and CH+ A4232 
[dashed line). Thick lines correspond to observed profiles while 
thin ones represent individual components. 



cause the blend with the stellar feature induces an underesti- 
mate of VFa(4232). We thus apply a small positive and con- 
stant offset to Wa(4232) measurements in order to get the same 
< N > for both transitions and finally compute the (a{N)^'^ 
weighted) average for each epoch. The best fit corresponds to 
SN/N ~ —0.65% per year and the 3(t bounds are —1.5 and 
0.27% per year (a non varying Af(CH+) is rejected at the 2a 
level); these values clearly exclude an increase as large as that 
seen for CH. 

Let us now compare the velocity profiles of CH and CH+ 
absorptions. Fig.|8]shows the apparent optical depth versus ve- 
locity for both species. It can be seen that no systematic shift is 
present; observed profiles are quite similar although the asym- 
metry of CH and CH+ lines is reversed. This is clearly seen 
when comparing the two narrow components needed to fit the 
profiles. (A comparison of the broad components would be 
meaningless because the CH+ one is ill-defined due to blending 
with the adjacent stellar absorption.) The CH and CH+ compo- 
nents have the same velocities but different CH+/CH ratios, as 
often observed (see e.g., [Crawford, 1989| and the ^ Per profiles 
given by Crane et al. 1995). 

3.3. Structure in CN 

Because CN is expected to respond strongly to density changes, 
the behavior of this species can give useful clues concerning 
the nature of the CH variations (Pan et al., 2001). The available 
observations were taken in 1993 (F93b), 1998 (H98) and 2002 
(F02). Unfortunately, when seen towards hot stars, the R(0) line 
is blended with stellar features from O II and C II (Fig.|9j see 
Meyer & Jura, 1985). As a result of the possibility of blending. 
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Fig. 9. The CN (0,0) lines in the F93b, H98 and F02 spectra. 
In the H98 spectrum, the S/N ratio is 700 per 0.07 A pixel. 
R(0) lines have been aligned at the rest wavelength. The three 
spectra have been normalised to 1 at A ~ 3875 and shifted in y 
for clarity. Note the broad stellar feature on the blue side of the 
R(0) Hnes. 



In the optically thin limit, the observed Wx ratio for R(0) 
and R(l) implies an excitation temperature T^xc = 3.6 ± 0.2 K 
(following the same method as Meyer & Jura, 1985) , signif- 
icantly above the value expected T^xc = Tcmb = 2.73 K if 
the excitation were due only to interaction with CMB radia- 
tion. As discussed above, the available data yield only loose 
constraints on 6(CN) and if the latter were small enough, satu- 
ration effects might affect the determination of Texc- To assess 
the importance of these effects, we compute Texc assuming b 
= 1, 0.5, 0.25, 0.20 km s^^ and get T^xc = 3.5, 3.4, 3.0 and 
2.8 K respectively. For lines of sight with similar extinction, 
6(CN) is observed to be about 1 km s^^ and values as low as 
0.2 km s^^ appear unlikely in comparison with available deter- 
minations (Crawford, 1995). Such low values for b would also 
be at odds with the results for CH A4300 ; no component in CH, 
that is associated with CN, has such a low 6-value. We conclude 
that Texc very likely exceeds Tcmb, implying that collisions 
with electrons contribute significantly in populating the N = 
1 level l IBlack & van Dishoeck, 199Tl l. If 6(CN) = 1 km s^^, 
Texc would be among the highest values quoted by Black & 
van Dishoeck (1991). 



we remeasured the line in a consistent manner In practice, this 
is not easy because the resolution differs from one spectrum to 
another and further, the broad stellar feature appears to vary in 
time. Results are given in Table lA.2l the larger value quoted by 
Herbig (1999) for the R(0) line is likely due to inclusion of the 
stellar feature. 

3.3.1 . Analysis of G. Herbig's spectrum 

Since the Keck data have the highest S/N ratio, we first discuss 
this spectrum and its implication on the CN column density, 
velocity distribution and excitation. The R(0), R(l) and P(l) 
lines are clearly detected, as shown on Fig.|9] 

To our knowledge, there exist no high resolution observa- 
tion with adequate S/N of the CN Unes towards HD 34078 that 
would allow us to measure their optical depth and investigate 
saturation effects in a direct way. Classically, the latter are esti- 
mated by comparing the relative strengths of the R(l) and P(l) 
lines (Meyer & Jura 1985). However, these CN absorptions are 
quite weak in the HD 34078 spectrum (W\ = 2.6 ± 0. 1 mA and 
M^A = 1 .5 ± 0. 1 mA respectively, while for the R(0) line, W\ = 
5.7 ± 0.2 mA), and the equivalent width ratio poorly constrains 
the b value. Indeed, only values smaller than 0. 1 km s^^ can be 
rejected (at the 2 a level); thus, within errors, the equivalent 
width ratio is consistent with CN lines being optically thin. 

Another way to estimate the optical depth of the 3874.6 A 
R(0) line is to compare its equivalent width to that of the R(0) 
line from the (1,0) vibrational band at 3579.96 A (Meyer et al., 
1989). In Herbig's spectrum, a weak feature is present at the ex- 
pected position with Wx = 0.55 ± 0. 15 mA; in the optically thin 
limit, a value Wx = 0.50 ± 0.02 mA would be expected. Again, 
within errors, the relative strengths of the two R(0) features are 
consistent with the assumption that these two lines are thin and 
only a lower limit on b can be obtained: b > 0.2 km s~^ (2 a 
limit). 



3.3.2. McDonald 2002 and 1993 spectra 

Very recently, the CN (0,0) lines have been reobserved at 
McDonald (R = 170 000) (Fig.|9j. The stellar feature around 
the R(0) line is not well defined due to limited signal to noise; 
apparently, it is weaker than in Herbig's spectrum. The R(0) 
and R(l) lines are clearly detected with Wx = 6.3 ± 0.5 mA 
and Wx = 2.6 ± 0.4 mA respectively. These values are com- 
patible with the 1998 ones. On the other hand, the 1993 (F93b) 
spectrum yields W^x = 5.1 ± 0.2 mA and Wx < 2.0 mA (2 
(7 limit) for the R(0) and R(l) Unes, implying Texc < 3.1 K. 
Altogether, the data suggest that both the iV = and iV = 1 CN 
column densities have increased, even if the variation is only 
marginally significant (AN/N = 12 ± 6%). Therefore, the CN 
variation could be as high as that for CH, but not much larger. 

3.4. Structure in DIBs Carriers 

The OHP spectra acquired for our study of CH and CH+ con- 
tain several strong DIBs. We now use these data together with 
older spectra to investigate the spatial distribution of DIBs car- 
riers. Some of them are known to behave differently with re- 
spect to the abundance of CH ( |Krelowski et al., 199 9). Then, 
an immediate goal is to determine whether the kind of small 
and large scale variations seen for CH are also present for DIBs. 

Among all detected DIBs, we initially selected those at 
5780, 5797, 6196, 6284 and 6614 A which are relatively strong 
(opacity reaching 6%) and narrow. For broader ones, the con- 
tinuum level and shape is difficult to determine and only large 
variations could be seen. Among the five features quoted above, 
we find that the sensitivity needed to detect variations with an 
amplitude comparable to that of the observed CH changes can 
be obtained only for the 5780 and 5797 A features. We then 
focus on these two DIBs. Since they are fully resolved (ex- 
cept possibly for some substructure), we can perform a direct 
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Fig. 10. Time variation of the 5780 A and 5797 A DIBs. Upper 
panels: the average of all OHP spectra (solid line) is compared 
to the spectrum obtained by Jenniskens et al. (1992) (dotted 
line). Lower panels: the ratio of the OHP and J92 spectra (solid 
line) is compared to the ratio obtained after artificially increas- 
ing by 10% and 15% the respective strengths of the 5780 and 
5797 A DIBs in the OHP spectrum (dashed line). 



comparison of the absorption profiles, once the spectra have 
been brought to a common wavelength scale. Equivalent width 
measurements would not be appropriate due to large continuum 
placement errors. 

We first discuss variations at the scale of 0.5 - 1 year (~ 
10 AU) among the spectra taken at OHP. Since the latter have 
been taken with the same instrument, the shape of the contin- 
uum adjacent to the DIBs considered is very stable, rendering 
the intercomparison much easier. We selected pairs of spectra 
(October 2000/March 2002 and February 2001/October 2001) 
in which Wa(4300) measurements show the largest difference 
(about 5.3% which coiTesponds to relative variations in iV(CH) 
of 7.2%). No significant variation is detected at the level of 10 
and 15% for the 5780 and 5797 A DIBs, respectively. To derive 
these upper limits, we artificially increase/decrease the strength 
of the DIB considered in the epoch 1 spectrum and then com- 



pare the latter to epoch 2 data. Then, the sensitivity reached is 
not large enough to detect variations as faint as those displayed 
by Af(CH) over a period of about one year. 

To study variations at larger scales, we use the spectrum ob- 
tained by Jenniskens et al. (1992) and compare it to the average 
of all our OHP spectra. The ratio of the two smoothed spectra 
is computed, and to ensure a good "relative normalisation", we 
fit values taken by this ratio in a few windows adjacent to each 
DIB with a polynomial of order 1 or 2. We then divide the first 
spectrum by this relative normalisation. The spectra thus ob- 
tained are shown in Fig.^Jand show no variation of the 5780 
and 5797 A DIBs at respective levels of 10 and 20%. Since both 
DIBs are optically thin, these upper limits are to be compared 
to the relative variation in iV(CH) which is estimated to be 20% 
over the same period. We then conclude that the 5780 A DIB 
at least has varied less than CH. 

3.5. Summary of observational results 

The main results that emerge from observations are the follow- 
ing: 

- over the past 12 years (or scales of 100-200 AU in the 
foreground cloud), the CH column density has increased by 
12-38%. A Hnear form (with dN/{Ndt) = 1.4% yr^^) or a step 
function ( AN/N = 22% in 1998) can both fit this overall varia- 
tion. Over the same period, the CH+ column density displays a 
markedly different behavior (a slow decrease with dN/ (Ndt) 
= -0.65%yr-i). 

- on the scale of a few months (~ 5 AU), correlated fluc- 
tuations in 7V(CH) and A^(CH"'") are strongly suggested by 
the data. A rough estimate of their amplitude over the pe- 
riod November 99 - March 02 can be inferred from the r.m.s. 
scatter of Wx measurements: (T(iV(CH))/iV(CH) w 3.6% and 
cr(7V(CH+))/iV(CH+) w 5.8%. 

- CH and CH+ velocity profiles show no significant time 
variation. 

- CN observations indicate a moderate increase of iV(CN) 
but a higher excitation temperature in 1998 and 2002, as com- 
pared to 1993. The recent T^xc value is significantly above that 
expected from excitation by CMB radiation alone. 

- The 5780 and 5797 A DIBs have been stable and the for- 
mer has varied by no more than half that seen in CH. 

3.6. Discussion 

3.6.1. AU scale structure and CH+ production 

We now discuss structure in the gas containing CH and CH"*" 
at the 1-10 AU scale as probed by the homogeneous observa- 
tions performed since 1999. We found indications for fluctua- 
tions of CH and CH+ absorptions: first, the scatter is relatively 
large and second, CH and CH+ fluctuations are apparently cor- 
related. Nevertheless, it should be stressed that these variations 
are quite small. Indeed, part of the observed scatter can be at- 
tributed to measurement eiTors and the real scatter in iV(CH) 
or iV(CH+) cannot exceed a few %. This is much smaller than 
the H2CO variations described by Moore & Marscher (1995), 
which attain 17% and smaller than N{\{ l) fluctuations seen in 
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the atomic phase over similar scales. In a simple model where 
most of the CH or CH+ gas is comprised of identical discrete 
entities, the number of such "clumps" must exceed 10^ to get 
fluctuations that weak. 

If the erratic variations seen since 1999 reflect the geometry 
of localised regions where CH+ molecules are produced, then 
it is natural to detect associated variations of CH since about 
30% of the CH molecules in a given sight line form from CH+ 
but the fraction approaches 100% for gas with densities less 
than about 100 cm"^ ( |Gredel et al., 2002| l. Presently, it is diffi- 
cult to check the consistency of this picture because the relative 
amplitudes of the CH and CH+ fluctuations are poorly con- 
strained by the available data. Alternatively, in a picture where 
CH"*" is produced throughout the low density medium (e.g., 
Draine & Katz, 1986, Pineau des Forets et al., 1986, Spaans, 
1995, Federman et al., 1996), the correlation between CH and 
CH"*" could be due to a physical mechanism which induces lo- 
cal density enhancements for all species. For instance, com- 
pression in regions where chaotic flows happen to converge has 
been invoked to account for the AU-scale structure in atomic 
gas (Jenkins & Tripp, 2001 [Hennebelle & Perault, 1999) and 
could also be effective in molecular gas. If the lifetime of these 
transient fragments is small enough, such a process would lead 
to comparable values for 5N/N(CH) and ^A^/A^(CH+); other- 
wise, the higher density might result in efficient CH"*" destruc- 
tion through reactions with H2 and rather produce an anticor- 
relation. One way to disentangle the above two scenarios in 
which small scale structure is either related to chemically ac- 
tive regions or just to pressure fluctuations would be to study 
the behavior of other species whose chemistry is not directly 
finked to that of CH and CH+. 

3.6.2. The long-term CH variation: evidence for a 
dense fragment ? 

Over larger scales (100 - 200 AU), CH and CH+ clearly dis- 
play distinct behaviors. Given the number of independent mea- 
surements and consistency among results, we consider that this 
difference is well established, even if the form of the CH in- 
crease is not well understood. In particular, we cannot exclude 
an abrupt variation by 20% occurring in 1997 - 1998. The cor- 
responding CN absorption has increased only modestly (6N/N 
= 12 ± 6%) between 1993 and 1998, but the data strongly sug- 
gest a change in excitation from T^xc < 3. 1 K in 1993 to T^xc = 
3.6 ± 0.17 K in 1998. The 1998 value is quite high in compar- 
ison with results quoted by Black & van Dishoeck (1991) for 
other lines of sight, implying a high electron density. The fact 
that the ratio A^(CN)/7V(CH) has remained roughly constant 
while A^(CH) and A^(CN) have increased, the ratio indicates 
that little change in gas density has accompanied the fluctua- 
tions in columns (see Pan et al., 2001). 

Consider a model in which the increase in A^(CH) by about 
20% over a time interval of 1 to 10 years is due to a newly in- 
tervening clump. A^(CH) is known to vary linearly with A^(H2) 
with a ratio 7V(CH)/iV(H2) « 4.0 lO^^ (Federman, 1982, 
Roueff, 2001), Towards HD 34078, a value larger by a factor 
of four is inferred (see Fig. 2 from Roueff, 2001). Assuming 



that such a ratio is still valid at the scales considered here, one 
can easily compute the associated increase in A^(H2): 1.1 10^° 
cm^^. Assuming an intermediate value of 50 AU for the trans- 
verse scale [A^(CH) has remained at its "high" value at least for 
4 years] and a comparable size along the line of sight, the im- 
plied density for that hypothetical fragment is around 1 .5 x 10^ 
cm"'^. In such a dense fragment, very little CH+ is expected to 
be present due to reactions with H2, accounting for the absence 
of an associated increase in A^(CH+). Similarly, the stability of 
the 5780 and 5797 A DIBs can be understood in this scenario 
since DIBs carriers are known to lie preferentially in diffuse 
media. However, the CN/CH ratio is expected to be large at 
these densities and the relative increase of A^(CN) should be 
much larger than that of CH, contrary to what is observed. The 
density quoted above could be lowered if the structure were 
elongated along the line of sight (cf. the model proposed by 
Heiles (1997) for atomic gas), but a very large (and unrealis- 
tic!) aspect ratio would be needed to bring it to values typical 
for that kind of material (a few times 10^ cm~^). In particu- 
lar, Federman et al. (1994) estimated ritot ~ 200 cm~^ from 
iV(CN)/7V(CH) and C2 observations, at a time when A^(CH) 
had its low value, but the column density ratio has not changed 
appreciably since then. While it is not clear which CH compo- 
nent contains CN as well, the cuiTent version of the chemical 
model (Knauth et al., 2001, Pan et al., 2001) gives similar den- 
sities. 

Another problem with this model (in addition to under- 
standing how the large overpressure inside the clump can be 
generated and maintained long enough) is related to the change 
in the excitation of CN suggested by the data. Since the vari- 
ation in A^(CN) has been quite limited (at most 20 to 30%), 
the higher Texc value estimated for the whole gas in 1998 and 
2002 implies an even larger T^xc in the clump if the rest of the 
gas has remained at T^xc ~ 2.7 K. Again, this points towards a 
large clump density or electron fraction. 

A last difficulty involves the observed velocity distribution. 
The presence of an additional fragment on the 2002 line of 
sight should have resulted in a new velocity component, but the 
whole distribution (in particular the two narrow components) 
has varied globally, without any noticeable change in shape. 
One could imagine that the additional component is very nar- 
row and does not show up at R = 200 000, but such a feature 
would be atypical with respect to profiles observed at very high 
resolution. 

3.6.3. Contribution from gas located near 
HD 34078 ? 

The difficulties encountered above when attempting to account 
for the CH - CN variations might be related to several pecu- 
liarities of the HD 34078 line of sight: i) the large amount of 
highly excited H2, ii) the large CH/H2 and CH+/H2 ratios (note 
however that the CH/CH+ ratio is in good agreement with the 
CH - CH+ correlation observed by Gredel (1997)), and iii) the 
large T^xc value for CN. Since point i) is attributed to the influ- 
ence of the UV bright star upon nearby material (Le Petit et al., 
2001, Boisse et al., 2002), one might wonder whether the prop- 
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erties observed for CH, CH+ and CN could also be the result of 
this interaction. In this model, CH, CH+ and CN are expected 
to be largely photodissociated in the PDR where H2 is excited 
to high energies. 

Then, to account for the observed CH and CH+ amounts, 
one may imagine that CH"*" production occurs in the PDR it- 
self, as has been suggested for the Pleiades (White, 1984) and 
IC 348 (Snow, 1993). In this scenario, CH+ is the result of 
reactions between C"*" and H2 in gas at elevated temperatures 
caused by enhanced H2 photodissociation. The large amount 
of CH then arises from the synthesis of CH+and the increase 
in Tpxc occurs because there is enhanced ionization produced 
by the star's UV radiation. One potential difficulty involves the 
velocity distribution which appears to be quite narrow (except 
for the broad component) and stable (cf. Herbig 1999). This 
suggests quiescent absorbing material rather than gas located 
near a bow shock around a 09 star moving at more than 100 
km s~^. Possibly, we are seeing "composite absorption" from 
both pieces of gas, one located close to HD 34078 and strongly 
influenced by it, and the other consisting of more distant mate- 
rial where the radiation field is enhanced by only a factor of a 
few relative to the average galactic field. 

3.6.4. Prospects 

Let us now suggest a few prospects concerning additional ob- 
servations that might clarify the situation. First, we wish to 
stress the importance of data other than those considered in 
this paper that could have been taken before 1999, especially 
spectra obtained in 1997 - 1998, which would be very useful 
in constraining the shape of the long-term variation in CH. We 
thus prompt any reader who would have such data in hand to 
contact us. Evidently, it is important to continue the CH and 
CH+ observations with increased S/N ratio in order to charac- 
terise better the properties of erratic fluctuations on timescales 
of about one year. At the same time, these data would give in- 
formation on the shape of the long-term variation of CH and 
CH+ absorptions. 

To date, FUSE observations have provided 4 spectra, all 
taken after 1999; their analysis indicates that A^(H2) variations 
- if any - are smaller than about 5% (Boisse et al., 2002). If 
A''(CH) keeps increasing, we shall be able to verify whether the 
CH - H2 correlation still holds at the 10-100 AU scales since 
the expected increase in A'^(H2) will exceed our detection Umit. 
The behavior of other species like C I and CO can be investi- 
gated also using these spectra. However, this can be done with 
better sensitivity with HST thanks to higher S/N and spectral 
resolution. HST/STIS spectra would provide invaluable infor- 
mation on the density via C I fine structure transitions and C2 
lines and allow us to follow with improved accuracy the small 
and large scale variations of many species (C I, H I, excited 
H2, CO, C2...). Repeating the C2 visible observations reported 
by Federman et al. (1994) could also give direct evidence for 
the high clump density quoted above. 
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Appendix A: Tables... 



Date 


Observatory 


S/N 


R" 


Transition'' 


References 


9-89 - >1-91 


KPNO" 


135 


85 000 


CH, CH+ 


Allen 94 


A91 


06-12/1991 


OHP (Aurelie) 


>100 


14 300 


CH, CH+(4232), DIBs 


JCllllisKClla Cl d.1., yZ, 




11-14-97 


1 Ci sivui 




ds nnn 

M-J \J\J\J 


CH CH* 


Krelowski et al., 99 


Ivy /a 


11 1 Q Q7 


ici sKUl 


Sf) Rf) 

JU ou 






Krelowski et al., 99 


Ivy / D 


ni 

yj L yj 






9nn nnn 








LZ,-yj 


iVlLJ_^UllalQ kjallUiiUlQ 


1 on 1 "^n 


6n nnn 

uu uuu 










iVlL.l-'UIlalQ CCllCllC 


iUU- i JU 


1 7n nnn 

1 / u uuu 






rUZ 


LU-yo 


JvcCK 


DUU 


Ac: nnn 
4j uuu 


prr prT+ pm 


Herbig, 99 




11 10 CiC\ 




1 Q ^ 1 /I n 
13J-14U 










12-22-99 




SI 0-^60 








R99b 


1-27-00 




120-200 








ROOa 


2-27-00 


OHP (Elodie) 


105-115 


32 000 


CH, CH+,DlBs 




ROOb 


3-18-00 




125-150 








ROOc 


10-01-00 




120-155 








ROOd 


2-10-01 




105-145 








ROla 


10-2-01 




105-145 








ROlb 


3-4-02 




150-160 








R02 



Table A.l. List of Observations. ° Resolution near CH A4300. Interstellar lines that have been included in this study. CH stands 
for CH A4300; CH+ for CH+ A3957 and CH+ A4232, unless specified; CN for CN A3875 (TablelO- DIBs stands for diffuse 
interstellar bands (see text for details). 





CH" 


CH-l- 


CH-^" 


CN" 


CN" 




4300 A 


3957 A 


4232 A 


R(0) 


R(l) 


A91 


48.3"" 


29.8^-" 


44.42-" 






J91 


42010.0 




52.0i"-°(b) 






r93a 


48.0^° 




44.01-^^ 






r93b 


52.01-2 


2930.8 


43.1"-^ 


5J"-2 


<2. (c) 


K97a 


49.6^° 


28.71-5 


39.8^-° 






K97b 


52.620 


25.42-" 


44.03-0 






H98 


58.0"-^ 


27.3°-2 


42 02 


5J"-2 


2.6"-i 


R99a 


55.7"-^ 


29.02-2 


38.02-0 






R99b 


55.4"-5 


27 91.4 


41.52-0 






ROOa 


55.2"-^ 


26.3I-5 


40.32-0 






ROOb 


57.81-2 


22 §2.8 


40.5^-° 






ROOc 


55.81-0 


2472.3 


40.02-° 






ROOd 


56.61-1 


23.21-'' 


39.02-0 






ROla 


58.8"-* 


31.92-^ 


44.03 






ROlb 


55.4"-^ 


24.41-3 


39.0^-" 






F02 


60.01-" 




44.02" 


6.3"-2 


2.6"--^ 


R02 


58.11-" 


25.31-^ 


42.02-0 







Table A.2. Equivalent width of interstellar lines (all errors are Iti error). ° Includes broad absorption (cannot recover the broad 
component). Does not include broad absorption from stellar line (visible in all spectra). 2a upper limit. 



